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INTRODUCTION 

This  paper  presents  a  method  for  predicting  the  reliability  (probability 
of  nonfailure)  of  basement  shelters  when  subjected  to  the  blast  effects  of  a 
single  nuclear  weapon  in  its  Mach  region.  The  method  is  described  with  refer¬ 
ence  to  a  reinforced  concrete  basement  shelter  whose  roof  slab  is  the  weakest 
structural  component.  This  is  generally  the  case  in  weak-walled  conventional 
buildings  when  the  first  floor  over  the  basement  is  at  grade  and  the  peripheral 
basement  walls  are  not  exposed  but  are  in  contact  with  the  soil.  In  such 
basements,  partial  or  total  collapse  of  the  slab  results  in  casualties.  Casu¬ 
alties  would  be  produced  by  debris  from  the  collapsed  slab,  the  building  above, 
and  by  pressure  build-up  within  when  the  shelter  envelope  is  breached.  The 
objective  then  is  to  determine  the  probability  of  roof  slab  collapse  and  on 
this  basis  to  determine  the  probability  of  people  survival.  The  paper  pre¬ 
sents  the  method  of  analysis  and  illustrates  its  application  by  means  of  an 
example  problem. 


STRUCTURAL  ANALYSIS 

The  form  of  structural  analysis  performed  is  described  in  Reference  (X). 
The  reinf^  ced  concrete  slab  is  modeled  as  a  single  degree  of  freedom  system 
whose  flexural  resistance  is  a  piecewise  linear  function.  The  resistance 
function,  see  Figure  1,  relates  the  flexural  slab  resistance  to  the  deflection 
at  its  midpoint.  Since  shear  is  a  possible  mode  of  failure,  the  analysis  is 
also  concerned  with  peak  dynamic  reactions  distributed  along  the  edges  of  the 
slab. 

The  blait  load  is  approximated  by  a  function  having  an  instantaneous  rise 
to  peak  overpressure,  followed  by  un  exponential  decay,  see  Figure  2.  It  has 
the  following  form  (Reference  2). 

F(t)  =  F] ( 1  -  t/td)e_t/td  (1) 

where  F^  =  peak  overpressure 

td  =  positive  phase  duration  or  the  overpressure 

The  spacial  distribution  of  the  blast  load  is  assumed  to  be  uniform  over  the 
surface  of  the  slab. 

Since  both  the  loading  and  resistance  are  complex  functions,  it  was 
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necessary  to  use  a  numerical  procedure  to  obtain  the  peak  midpoint  deflection 
and  the  peak  dynamic  reactions.  The  primary  equations  used  in  the  analysis 
are  the  following. 

KLmV  +  R(y)  =  F(t)  (2) 

where  =  the  load-mass  factor  (Reference  ]_) 

M  =  the  total  mass  of  the  slab 

R(y)  =  flexural  resistance 

F(t)  =  load=time  history,  se*  Equation  (1) 

V(t)  =  C1R(y)  +  C2F(t)  (3) 

where  V ( t )  =  the  dynamic  reaction  along  the  given  edge  (a  or  b)  of  the  slab, 
see  Figure  3 

C^,C2  =  constants  whose  values  depend  on  the  aspect  ratio  of  the  slab 

PROBABILITY  OF  FAILURE 

In  the  case  of  two  failure  modes,  the  probability  of  slab  failure,  P(F), 
is  (Reference  3) 

P(F)  -  1  -  [1  -  P(Fb)][l  -  P(Fy)]  (4) 

when  the  modes  are  independent,  and 

P(F)  =  max  [P(Fb),  P(Fy)]  (5) 

when  the  modes  are  highly  correlated.  In  Eqs.  (4)  and  (5),  P(F.  )  is  the 
probability  of  failure  due  to  flexure,  and  P(F  )  is  the  probability  of  failure 
due  to  shear.  The  actual  probability  of  failure  is  between  these  two  bounds. 

Probabilities  of  failure  due  to  flexure  and  shear  were  each  computed 
using  the  following  expression 


P(F)  =  1  -  H - MLirZl) - }  =  i  .  $[ML]  (6) 

/  inLci  -  rspri  - 

where  4>(  )  =  the  cummulative  density  function  of  the  standard  normal  distri¬ 
bution 

r  =  the  median  value  of  the  resistance  parameter  in  flexure  or  shear 

s  =  the  median  value  of  the  load  parameter  in  flexure  or  shear 

fir  =  coefficients  of  variation  of  the  resistance  and  load  parameters 
respectively 
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0  =  median  safety  factor 


=  total  degree  of  dispersion  of  the  safety  factor 


For  the  case  of  flexural  response  the  median  safety  factor  u  is  taken  as 
the  ratio  of  y  /y  ,  where  y  is  the  ultimate  (collapse)  midpoint  deflection 

Jm  Jp  ■'m  r  I  r 

of  the  slab,  see  Figure  1,  and  yp  is  the  midpoint  deflection  at  a  given  load. 
The  value  of  y  is  taken  as  0.15a,  where  "a"  is  the  short  span  dimension  of 
the  slab  (Reference  4) . 


For  the  case  of  shear  response,  the  median  safety  factor  is  taken  as 

the  ratio  of  v  /v  ,  where  v  is  the  ultimate  unit  shear  capacity  of  the  slab 
m  p  m 

and  Vp  is  the  corresponding  shear  stress  at  a  given  load.  The  shear  stress 

is  computed  at  the  periphery  of  the  slab  by  the  use  of  dynamic  reactions 
mentioned  earlier.  The  ultimate  unit  shear  capacity  of  the  slab  is  based  on 
the  following  formula  which  is  the  standard  ACI  (Reference  5)  formula  modi¬ 
fied  as  suggested  in  Reference  6. 


v  =  1.5(2  f'  ) 
m  dc 


(7) 


where  f^  =  1.25fl  =  the  ultimate  compressive  strength  of  concrete  increased 

to  account  for  the  increase  in  strength  due  to  dynamic  loading 
conditions  (Reference  1_) 

SAMPLE  APPLICATION 

Figure  3  shows  the  plan  view  of  a  reinforced  concrete  slab  whose  rein¬ 
forcing  steel  extends  over  and  beyond  the  supports.  Supports  are  continuous 
along  the  edges  of  the  slab.  The  reinforcement  in  the  short  direction  is 

0 . 27 ( i n ) 2/f t  (572mm2/m)  and  in  the  long  direction  is  0.19(in)2/ft  (402mm2/m). 
The  slab  is  9-in  (228.6mm)  thick.  The  compressive  strength  of  concrete,  f1 
=  3000  psi  (20.7  MPa)  and  the  yield  strength  of  reinforcing  steel,  f  =  c 
60,000  psi  (414  MPa).  y 

In  performing  the  analysis,  the  following  parameters  were  treated  as 
random  variables,  i.e.,  F-| ,  t^,  f^,  f  ,  A$  (cross-sectional  area  of  rein¬ 
forcing  rods),  d(effective  depth  of  the  slab). 

Coefficients  of  variation  of  the  basic  parameters  were  obtained  from 
available  experimental  data  (Ref.  4,  ]_,  8).  Corresponding  coefficients  of 
variation  of  slab  resistance,  peak  deflection  and  peak  shear  stress  were 
determined  on  the  basis  of  a  first  order  approximation.  Reference  3. 

This  slab  was  analyzed  when  subjected  to  a  series  of  blast  loads  of  in¬ 
creasing  intensity  with  durations  corresponding  to  a  1-MT  surface  burst. 
Results  of  the  analysis  are  shown  in  Figure  4-  and  Figure  5.  Figure  4  shows 
the  probabilities  of  failure  in  flexure  and  shear  taken  separately  and 
determined  on  the  basis  of  Eq.  (6).  Figure  5  shows  the  bounds  on  the  pro¬ 
bability  of  failure  computed  on  the  basis  of  (4)  and  (5). 
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SUMMARY  AND  CONCLUSIONS 


A  method  for  predicting  the  probability  of  failure  of  structures  by  con¬ 
sidering  multiple  failure  modes  was  formulated.  It  was  applied  to  the  analy¬ 
sis  of  a  reinforced  concrete  slab  when  subjected  to  a  uniformly  distributed 
blast  load  over  its  surface.  Currently  available  criteria  for  failure  due 
to  flexure  and  shear  (Ref.  _4  and  6)  were  used  in  predicting  the  probability 
of  failure. 

This  method  is  capable  of  considering  all  major  components  of  a  structure, 
the  respective  failure  modes  of  each  component,  and  of  predicting  the  pro¬ 
bability  of  failure  of  the  structure  as  a  whole. 
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